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Abstract 

This paper presented the existence and representation of solution to the 
■ initial value problem for the general n-term linear fractional differential equa- 

\ tion with generalized Riemann-Liouville fractional derivatives and constant 

C*") ' coefficients by using operational calculus of Mikusinski's type. We proved 

that the solution of the initial value problem exists if and only if some initial 
2^ ' values should be zero. 
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1 Introduction 

In recent years there have been considerable interests in the theory and applications 
of fractional differential equations (JT] [3] |TT] [T2] [14] [T2l )- Indeed, it is well known 
that fractional derivative is excellent tool for description of memory and heredity 
effects ifTUll . It is therefore important to understand the advantages of using frac- 
tional derivative in classical equations. 

Applications of fractional calculus require fractional derivatives of different 
kinds O. There exist several different definitions of fractional differentiation 
0Q21. The (right-sided) Riemann-Liouville fractional derivative is defined by 

D o + m ■■= ^srvo*) > * > o (i) 
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where 



IoV(z) := J-r /% - tT^f^dt , x > (2) 

ig + /(z) := /(z) , x > 

is the (right-sided) Riemann-Liouville fractional integral of order a with lower 
limit 0(see e.g. CJOl). The (right-sided) Caputo fractional derivative is defined 
by 

d n 

c D« + f(x):=lZ-<*—f(x),x>0 (3) 

whenever the right-hand side exists(see e.g. fTlfTBl). 

In [4] the generalized Riemann-Liouville fractional derivative (GRLFD) of or- 
der a and type (i is defined as 

Dgfix) := (lSt a) ^+ P)in ~ a) f))^) > W 

where n — 1 < a < n G N and < /3 < 1, whenever the right-hand side exists(see 
e.g. (H). This generalization gives the classical Riemann-Liouville fractional dif- 
ferential operator (1) if /3 = 0. For /3 = 1 it gives the Caputo fractional differential 
operator (3). The concept of GRLFD is appeared in the theoretical modeling of 
broadband dielectric relaxation spectroscopy for glasses 0. Some properties and 
applications of the GRLFD are given in ffl Q2] [13 OH 23 . Several authors 
(see |2j [151 [191) called (4) the Hilfer fractional derivative or composite fractional 
derivative operator. 

More recently, new results for the fractional differential equations with GRLFDs 
were obtained in |f2] [15] H2J [20]. In [20] the authors provided an approach based 
on the equivalence of a nonlinear Cauchy type problem with GRLFD to a nonlin- 
ear Volterra integral equation of the second kind in spaces of summable functions 
on a finite interval of the real axis and proved uniqueness and existence of the 
solution using a variant of the Banachs fixed point theorem. In [19 ] they investi- 
gated the solution of spacetime fractional diffusion equations with a generalized 
RiemannLiouville time fractional derivative. 

The Mikusinskis operational calculus was successfully applied in ordinary dif- 
ferential equations, integral equations, partial differential equations and in the the- 
ory of special functions. For example, in JS] this operational calculus has been 
applied to solve a Cauchy boundary-value problem for a certain linear equation 
involving the Riemann-Liouville fractional derivatives. In [9] by using the opera- 
tional calculus of Mikusinskis type for the Caputo fractional differential operator 
they obtained exact solutions of an initial value problem for linear fractional dif- 
ferential equations with constant coefficients and Caputo fractional derivatives. 

In J6l an operational calculus of the Mikusinskis type was introduced for the 
generalized Riemann-Liouville fractional derivative operator and it was applied to 



2 



solve the corresponding initial value problem for the general n-term linear frac- 
tional differential equation with constant coefficients and GRLFD of arbitrary or- 
ders and types. But, in [6] they provided initial values independent on the type of 
generalized Riemann-Liouville fractional derivatives in n-terms linear fractional 
differential equation(see equation (37) of [6]) and defined the space of solutions 
independent on the type(see definition 1 of (H). They did not provide a necessary 
and sufficient condition for existence of the solution of n-terms linear fractional 
differential equation and thus they made some mistakes^. For example, accord- 
ing to the formula of solution given in ||6l (see equation (42) of [6]), the solution 
of the following initial value problem 

D™>°y(x) - D° fy(x) = 0, l ° + 3 y(0) =C„eR, / ° + 5 y(0) = C 1 € R 

is y(x) = (Co + C\)x ' 2 Eq.2,o.7(% ' 2 ) » but this is not a solution of the initial value 
problem. According to our result, the solution exists if and only if the initial value 
C\ should be zero. 

In this paper, we present the existence and representation of solution for the 
n-term linear initial value problem with generalized Riemann-Liouville fractional 
derivatives and constant coefficients by using operational calculus. We start with 
some preliminaries and operational calculus for the generalized fractional deriva- 
tive in Section 2. In Section 3, we provide a necessary and sufficient condition for 
existence of the solution of the initial value problem and some examples. 

2 Preliminaries and Operational calculus for generalized 
fractional derivatives 

Definition 1. A function f(x), x > 0, is said to be in the space C™, m € No = 
N U {0}, if there exists a real number , such that 

f(m) {x)=x p fi{x) (5) 

with a function fi(x) in C[0, oo). Especially the space is denoted by C a . 

Clearly, C a is a vector space and the set of spaces C a is ordered by inclusion 
according to 

C p C C a O a > p. (6) 

Lemma 1. /[9]/ The Riemann-Liouville fractional integral 1 'q + , a > 0, is a linear 
map of the space C 7 , 7 > — 1, in to itself, that is, 

-^0+ : Cry — > C a + y C Cj. (7) 

It is well known, that the operator Iq + , a > has the following convolutional 
representation in the space C 7 , 7 > —1. 

/ +/(x) = h a {x) * f(x) = / f(x)h a (x - t)dt , x > 0, (8) 

J 
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x a-l 



K{x) := — - , x > 0, a > 0, (9) 
r(a) 

where x 

(g*f)(x):= [ g(x-t)f(t)dt, x>0, (10) 



is the Laplace convolution. Moreover, the following properties of the Riemann- 
Liouville fractional integral are well known 

w& /(*) = iS+i&n*) = C 3 . /3 > o ; 

(/ a + .../ V)(x) = (/ T/)(^)- (11) 

S v ' 

n 

Lemma 2. ^ Lef / iraC^, m G N , /(0) = ••• = /(^(O) = and 
5 € C?£]_. 77je« the Laplace convolution 

h(x) := [ X f(t)g(x-t)dt (12) 

JO 

is in the space C m x and h(0) = ■■■ = /i( m )(0) = 0. 

Now we defined the following space of functions and use in operational calcu- 
lus. 

Definition 2. Afunction y G C_x is said to be in the space £1% ifl^. /3 ^™ 
Cl-y, n- 1< a < n G N/ora > 0, Qleq/3 < 1. 

Then the following lemmas hold. 

Lemma 3. Lef y G f2^, re — 1 < a < n G N. /ze?i f/ie Riemann-Liouville 
fractional integral (2) and the generalized fractional derivative (4) are connected 
by the relation 

(I$ + D«fy)(x) = y(x) - y a #{x), x > 0, (13) 

where 



k=0 



^— ' r(fe — n + a — pa + pn + 1) ax' ^ 

(14) 



As in the case of Mikusinskis type operational calculus for the Riemann-Liouville 
or for the Caputo fractional derivatives we have the following lemmas (see e.g. 
HIE)) 

Lemma 4. The space C_i with the operations of the Laplace convolution and 
ordinary addition becomes a commutative ring (C_i, *, +) without divisors of 
zero. 
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This ring can be extended to the field M_i of convolution quotients by follow- 
ing the lines of the classical Mikusinski operational calculus; 



M_! := (C_i x (C_i \ {0}))/ ~, (15) 
where the relation ( ) is defined, as usual by 

(/, g) ~ (A, 91) & (f*gi)(x) = (g*fi)(x). 

This relation is an equivalence relation in terms of Lemma 4. For the sake of con- 
venience, the elements of the field M_i can be formally considered as convolution 
quotients The operations of addition and multiplication are the defined in M_i 
as usual: 

/ . h f*gi+g*h n „ 

— I := (lo) 

5 gi g*gi 

and 

L . ll ■= f * f 1 (17) 
5 91 ' g * 91 ' 

Then the space M_i with the operations of addition (16) and multiplication 
(17) becomes a commutative field (M_i, -, +), where the unit element is de- 
fined by I := ^,and the zero element is defined by := j, / / 0. The flied 
(M_i, •, +) is the quotient flied of the ring (C_i, +). The ring (C_i can be 
embedded into the field M_i by the map: 

/ (/, J) = a > 0. 

= t-2-, a > 0. (18) 



h a h a starh a h 2a 

Lemma 5. For the function y € C_i, Riemann-Liouville fractional integral 
of y can be represented multiplication in the field M_i : 

(ig+y) (*) = -■*/. (19) 

s a 

A«<i generalized fractional derivative D^y of the function y G O?, a > 
0, < /3 < 1, ca« &e represented as multiplication in the field M_i of convolution 
quotients 

(D^y)(x) = s a -y- s a -y at p. (20) 
Formula (11) means that for a > 0, n G N 

/i™(x) := underbraceha ★ ■ ■ ■ ★ /i Qn = h na (x). 
This relation can be extended to an arbitrary positive real power exponent: 

:= h Aa (x). 
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For any A > 0, then inclusion h^(x) € C_i hold true and the following relations 
can be easily proved 

hP ct -kh? a = h/3 a -k h 7a = h^ +1 ) a = /i^ +7 , bP ai = W a% <^> azi/3 = a 2 j. (21) 

The above relations motivate the following of the element s a with an arbitrary real 
power exponent A. 



h- x ifA < 



I ifA = (22) 
It if A > 



For any a, a%, a 2 , /3, 7 > 0, it follows from this definition and the relations (21) 
that 

*a ■ S Z = S(/3+ 7 )a = S a +7 > S «i = S l 2 ° "l/ 3 = "27- 

For the application of the operational calculus to solution of differential equations 
with generalized fractional derivatives it is important to identify the elements of the 
field M_i, which can be represented by the elements of the ring C_i. One useful 
class of such representations is given by the following lemma (see e.g. lU|9l). 

Lemma 6. Let /3 > 0, a-i > 0, i = 1, ■ ■ • ,n. Then 

l£i — r P a i Ei \ a (Ai r aia ■■■ A T a " a ) (2Vi 

with the Mittag-Leffler function 

E( ai ,...q n ),b(zi, ■ ■ ■ , z n ) := (Mi,--- '^Ylb + Y™ ' a-l)' 

k=oh+-+i„=kh,-,in>o ^ L>j=i i i> 

aj > 0, 6 > 0, Zj E C (24) 

multinomial coefficients 

k\ 

{k, l\, • • • , Z n ) := 



M x ••• x L! 



When a € Z and 6 £ Z satisfy a < b, Z 0> & is used for the set of all integers i 
satisfying a < i < b. 



3 Existence and uniqueness of the solution 

In this section the constructed operational calculus is applied to solve linear frac- 
tional differential equations with generalized derivatives and constant coefficients. 
We consider the linear differential equation 

n 

D a Q f°y{x) - Y, a t D^y(x) = g(x) (25) 
i=i 
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and the initial conditions 

_ / (l-ft)(m i -a i ) y(0+) = yk _ {l _ M{mi _ ai) G R; i G Zon) fc G Zomi _ 1; 

(26) 

where a.; > 0, m; L — 1 < a« < m 8 G N, < < 1, i £ Zo, n and ^ £ R, i G 
Zi, n. The ordering 

m - (1 - /3 )(m - a ) > ■ ■ ■ > m„ - (1 - (3 n )(m n - a n ) (27) 

is assumed without lost of generality. 

Definition 3. The function HILo ^B- ™ ca ^ ec ^ tne solution of the initial value 
problem (25) and (26), if y satisfied the equation (25) and the initial value (26). 

Theorem 1. The initial value problem (25) and (26) with homogeneous initial 
conditions 

^_ / (i_-ft)(^-«i) J/(0+) = o, i G z 0in , fc G Z , mi _i, (28) 

where g(x) G | ^J 1 7^° ^ iv has a solution y(x), unique in the space {y^ =Q Q^ 1 , 
of the form 

y ( x )= /V- 1 S (ao _ ai) ... ao _ an)jQO (ait a «'- ai .--- ,a n r°- Q ") 5 (x-i)^ (29) 

JO 

Proof: Let G P|?=o ^ft and satisr i ed (25) and (28). Then the following 
algebraic equation in the field M_i of convolution quotients is obtained 

n 

1=1 

from homogeneous initial value problem (25) and(28) by using the Lemma 5. This 
linear equation has a unique solution in the field M_i 



y = — —^n 9 (31) 

2-a=l a i s cti 



From (23) we get 



i=l a i s cti 1 — Z-a=l a i s -(a -ai) 
X E(ao—ai,---ao—a„),a () ( a 'l x ° S''' -, a n% ° ) (32) 

Therefore we obtain 

= x a °" 1 £; (ao _ ai) ... ao _ Qt>))QO (aix a '>- ai , ■ • • , a n x a °- a ") *g(x). 
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Thus, we have the solution (29). Now we consider y{x) € HiLo ^ft • Let y lk := 
Io+V for 7fc = (1 - Afc)( m fc - a k ), k € Z ,„. Then we have 

= h lk {x) * (x Q °- 1 £; (Q0 _ Ql) ... Q0 _ aB)ia() (aix Q0 - ai , ■ • • ,a n x a ^))*g(x) 

= (x ^(Q -ai,-ao-C(„),ao (o>lX , <2 n X ")) -k (hj k (x) -k g(x)) 

= / t E( a() _ ai,---ao~a n ),ao ( a l^ ° 1 j ' ' ' i a nt ° n )9~{k( X ~ t)dt 
JO 

where g lk {x) := l2+g{x) G | ^J 1 ^ • Whereas, from (31) we have 

y(x) = i£ 5 (z) + «^o 7 r 7l y(^) ( 33 ) 
i=i 

and 

n 

i=i 

By (34) we obtain 

y lk {x) = (/jr mma Vi)(x), ^ € { ^ ^ JJ (35) 

Combining now the relations (34) and (35) and repeating the same arguments p 
times (p = [ao/(ao — a{)] + 1), we arrive at the representation 

2/7*0*0 = Uo+V'p)^), YV € | ^ij ^ (36) 

In the case ao = m o £ N, it follows from (36) that 

y<7> = ip p £ d, y 7o (0) = • • • = = °- 

If ao ^ N then since 

M*) = jL_ e c*r\ ^ao(o) = • • • = ^?- 2) (o) = o. 

we obtain by using Lemma 2. In the case similarly with above argument, we have 

y, k eC^,y, k (0) = --- = y^- 1 \0) = 0. 

Hence is a unique solution of homogeneous initial value problem (25) and (28). 
The proof of theorem is completed. 
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Theorem 2. For the initial value problem (25) and (26) with inhomogeneous 
initial conditions, in the case (1 — /?o)(^o — a o) = • • • = (1 — /3 n )(m n — On) = 7» 
?/ie?i f/ie initial value solution(25) and (26) has the unique solution 

y(x) = y 9 {x) + ES 1 Vk-^Tif^T) + 

n 

4_ \ n .rpk-J+ao-ai p ( a -ai ... a -a n u 

' / j u f h ^(ao— ai,---ao— a n ),k+l — 7+ao— «i ^ "l- 1, > > "n-^- 

(37) 

m f&e space HiLo ^ft' wnere h & No, fe = 0, ■ ■ ■ , mo — 1 are determined from 
the condition 

m h> k + l (38) 
m h+1 <k 

in the case rrii < k, i = 1, • • ■ , n we sef := 0, and if mi > k + 1, i = 0, ■ ■ ■ , n, 
tfjen / fc := n. 7/7 = (1 - /3 )(m - a ) = ■ ■ ■ = (1 - A - l)(mj - 1 - a/ - 1) ^ 
(1 — [3i)(mi — a/), < / < n, then there exists a solution y(x) in the space 

n™=o ^ft> tf and on fy if 

yk-(l-/3i)(mi-ai) = 0, i G Z/ jn , A; G Z , mi _i, 

y fc _(i_ ft ) (mi _ Qi ) =0, i £ Z 0j /_i, fc G Z ,m, M = [m / -(l-/3 / )(m / -a / )-l+7], 

(39) 

and a unique solution y[x) is represented in the form . 
y(x) = y 9 {x) + 2~2T=m yk-A r(k-7+i) 

n 

_i_ \ „.J-7+«o-«. p, , , _ ( ao-ai . .. T «o-««\l 

t / j u, « x ■ L -'(ao—ai,---ao—a„),k+l—j+ao—ai\ u 'l Jj > ) "n-^ // 

(40) 

where G No, fc = 0, • • ■ , mo — 1 are determined from the condition 

m h > k + 1 .... 

(41) 

mi k +i < k 

in the case mi < k, i = !,■■• ,nwe set := 0, and if mi > k + 1, i = 0, • • • , n, 
?/ie« /fc := n. //ere z's a solution (29) of homogeneous initial value problem 

(25) and (28). 

Proof: By using lemma the following algebraic equation in the field M_i of 
convolution quotients 

n n 
S a ■ y ~^2 a i s a t ■ y = 9 + S ao ' ^aoA ~ ^2 ° iSa » ' Va ^ ^ 
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from inhomogeneous initial value problem (25) and (26), where y a% fi l is repre- 
sented by (14). This linear equation has a unique solution in the field 

II n 

V = - J™ " 9 + - 7™ — ( s a ■ Vao,l3 — ^2 aiS ai ■ y ai ,Pi) 

(43) 

Here y a := ^ • q is represented by (29) and is the solution of homo- 

y s otQ— Z-,i=l a i s ai 

geneous initial value problem (25) and (28). 
Let 



Vh ■= — ™ (s ao ■ Vao^o ~ ^ a * S "« ' V<*i,Pi) ( 44 ) 

s ao 2^i=l a i s ai i=1 

we have 

Vh ~ ^-Er=i a i*«i-«o ^ k=0 yk-{l-Po)(m -a )Z-(k-(l-l3o)(mo-ao)+l)- 
n mo — 1 

-^aiS^-ao ^ y fc -(l- / 3 )(mo"ao) ,S '-(fc-(l-/3o)(mo-ao)+l)} (45) 



i=l fc=0 

In the case (1 — Po)(mo — ceo) = • • • = (1 — f3 n )(m n — a n ) = 7 we have 

mo — 1 r T sr^h mo— 1 t 

E J J — Z^i=l a i s aj-ao rri ^ 
t -sr n = 1^ Vk -t { I+ TZ. V" 



mo-1 fc _ 7 n 

= E Vk ^ T (k-y+l) + ^ «^"" 7+Q °" a ^ a ^^ + l-7 + ao-a^(x)} 
fc=0 ^ ' j=/ fc +i 

from (45), where & = 0, • • • ,mo — 1 are determined from condition (38). 
According to the definition of the numbers we have m; < fc for z G Z; fe+l n . It 
follows then, that A; — 7 + ao _ > a o — gamma, i G lk + 1, n. Using this 
inequality, we readily get the inclusion Iq + G C m " and y(x) G n™ =0 f^\ For 

n fc := r^f^ry + Ztl k+ i a^^ 7+a °- a ^ a ^(.),fc+i-7+«o-a^(^) we have 

(W (0 (0)={ J j^jj ,fc,I = 0,- > m o -l. 

Therefore satisfies the initial condition (26) and hence y(x) = y g (x) + yft(x) is 
a unique solution of the initial value problem (25) and (26). If 7 = (1 — A)) (mo — 

oo) = • • • = (1 - A - l)(mj - 1 - ai - 1) ^ (1 - - aj). < Z < n, then 

we have 



mo — 1 I— 1 m;- 1 

= { E ^-7 S -(fc-7+l) — E a * S "i-ao E yfc-7 S -(fc-7+l) — 

n m;- 1 

-y^ajS a ._ Q0 ^ yfc- 7 -(l-ft)(m i -a i ) ,S -(fc-7-(l-ft)(nii-ai)+l)J' = 



i=Z fc=0 



10 



m -l j _ v^'fc 

_ 1 2^=1 a i s ai -a<) 

~ Z> yfe-7 S -(fe-7+l) j_ yn 

k=0 l^i=l 1 a 



ao 
rrii — l 



j _ ST^ri — '^ / a i s a i -a E yfc-7-(l-/3 i )(m i -a i )' s -(fc-7-(l-ft)(mi-Q i )+l) 

1 ^i=l°< s ai-ao i= , fc=0 



from (45), where 4 is defined by (41). From above equation we have 

mo — 1 n 

7 n V I ' ) i 1 



k=0 ' ' ' i=l k + l 

n nii — 1 



E fli E Vk-(l-(Wm i -a i )X ao ° l+k (1 ft)(m ' ai) J E;(.) 1 ao-a i +fc+l-7-(l-|9 i )(m i -a i )( a: ) 



i=/ fc=0 



In order that y^(x) G HiLo > tne following relation should be satisfies 

a - on + k + 1 - 7 - (1 - A) (m, - a*) + 7 - m = 
= fc + 1 - m, + /3j(mj - «i) - /3 (m - a ) > 

for every i G Zj n , A; € Zo, mj -i. This relation satisfies only in the case A; = rrii — 1 
and ft(mj - a*) > /3 (m - a ). Therefore y jfc _ (1 _ /3 . )(m ._ Q . ) = for every 
i G Z/ ra , At G Zo imi -2 and, as yh{x) must satisfies the initial condition (26), we 
get j/mi-i-fi-ftx^-oi) = for every i G Zj )fl . Also, as the relation 
A - 7 + (1 - /3/)(m; - a/) - fa > 0, fc = 0, • • ■ , m - 1, fa = 0, ■ ■ ■ ,m/ - 1 
must satisfies , the following relations should be satisfies i/fc_ 7 = for A: = 
0, • • • , M, where M = [mj - (1 - A)(™z - ai) - 1 + 7]- 
Therefore we obtain 

mo— 1 fc— 7 n 

y fc (x)= E ^( fffc^Tl) + S «^ fe ~ 7+a °~ a ^ a ^(-), fc+ l-7 + «o-a^(x)) 

k=M \ ™ ' i=; fc +i 

Hence the solution is y(x) = y 9 (x) + yh(x). The proof of theorem is completed. 

Example 1. If (3q = ■ ■ ■ = /3 n = 1 then the initial value problem (25) and (26) 
becomes the initial value problem with Caputo fractional derivatives 

n 

c DZv{x) - E <D^y(x) = g(x) (46) 
i=i 

d k 

^y(0+) = y k G R, k G Z , mi _i, (47) 

where ao > a± > ■ ■ ■ > a n > 0, rrii — 1 < on < mj G N, i G Zo, n > ^ n <^ ^ e 
space PlILo ^ coincides with C™°. From theorem 2 the initial value problem 
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(46) and (47) has a unique solution 

rx m 0-l k 

n 

_i_ \ k+OLo-ctiT? ( ao-ai . .. ao-a„\\ 

(48) 

where k = 0, • • • , mo — 1 are determined from condition (38). 

Example 2. If /3q = ■ ■ ■ = (3 n = then the initial value problem (25) and (26) 
becomes the initial value problem with Riemann-Liouville fractional derivatives 

n 

D^y(x) - J2 <HD$.y{x) = g(x) (49) 

i=l 

d k 

-p;I™+ ai 2/(0+) = y k -( mi - ai ) Z ,„, k G Z , mi _i, (50) 

where ao > ct\ > • • • > a n > 0, m,j — 1 < a% < rrii G N, i G Zo, ra • 
mo — ao = • • • = "^n — «n = 7 ^en f/ie initial value problem (49) and (50) has a 
unique solution 



y(x) 



r^°- 1 ^(ao-a 1 ,..ao-a„),ao(aii a0 " Q1 >- " " ,a n t a ^)g(x - t)dt+ 

Jo 



mo-1 ^ fc _ 7 



fc=o v / 1 / i=i fc +l 

(51) 

(/i<?re £ = E {ao _ au ... ao _ an)jk+1 _ 1+aQ _ a . (aix a °~ ai , • • • , a ra x a °- a ™)) in the space 
fXi=o ^o* = ^o'- ^ ere 'fc> A; = 0, ■ ■ ■ , mo — 1 are determined from condition (41 ). 
If 7 = mo — «o = • • • = mi-i — a^_i ^ mi — ai, < I < n then exists a solution 
y(x) in the space HILo ^o' = ^o*> if an d only if 

Vk-imi-ai) = 0, i G Zj jTl , A; G Z , mi -l 

y fc _( mj _ Qi ) = 0, i G Z ,j_i, fc G Z ,m, M = [mj - (m z - a*) - 1 + 7] (52) 

Then, the unique solution of the initial value problem (49) and (50) is represented 
in the form of 

V(x) = r^-^o-m.-ao-a^oM 00 "" 1 '--- ,a n t a °- a ")g(x-t)dt+ 
Jo 

mo— 1 fc— 7 n 

k=M V ' 7 i=« fc +l 

(53) 
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Here E=E {aQ _ au ... ao _ an ^ k+1 _ 1+ao _ cti (a 1 x a ° ai ,---,a n x a ° Q ") and l k , k = 
0, • • • , mo — 1 are determined from condition (41 ). 

Example 3. We consider the so-called composite fractional relaxation equation 
(see [JU) 

rxj t f{t) + r?D%ff(t) + f(t) = (54) 

/(0+) = 1, (4^ 1 - Q )/)(0+) = (55) 

where < a < 1, < /i < 1, < n, < oo. By using theorem 2 the initial value 
problem (54) and (55) have a unique solution f(t) in the space Q,q f] $7^ f] $7q, 
and only if 

(4^ )(1 ~ Q) /)(0+) = = 0, (56) 

Then the unique solution is 

f(t) = £ (1 _ aj 1)f U-^t 1 " , - -^—t), t > 0. (57) 
n r — 1 
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